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ABSTRACT 

The premium calculation in insurance is of great 

importance. Knowing past claims L , .... | one wishes 

1 n 

to forecast the next claim £ , , or in a more general 

fashion, forecast a function of £ , f(| ) , 

n+i 'n+1 

Since we are interested in a least squares forecast of 

f(.| ) , the problem readily reduces to calculating 

E[f(£ , , ) I C, , .--j £ J • This can formally be done by use 
n+1 1 n . 

of Bayes rule. As Bayes rule is difficult to apply the 
theory of approximation with an orthonormal set of functions 
is exploited in this paper, to reveal the nature o£ our 
problem. The credibility premium and credible distribution 
follow as partial results after a generalization of Ericson's 
Theorem. 
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INTRODUCTION 

Let C_ = £-, s -.-s i i.i.d. random variables with density function 
P 5 | e (x | 6) , i.e., C ± ~ P c | e (x | 9) • 

Let their joint density function be p i (x | 0) , Since they are independent 
n 

p |e (x 1 e) = n P |Q ( X . | 0) . 

Also let (where can be a vector) be distributed as p Q (8) . We 



will call this density the prior. 

Now suppose that we have observed E, , . . . , E, (which correspond to the 

individual experience in insurance terminology) and we want to find a function 

<j>(jj) to approximate in the least square sense a function of the next observation 

f(£ ,,) . The above is summarized as: 
n+i 

Find (f)(i) so that 



I = E 



r J f (5 n ,-,) - + <£)] is minimized, 



'n+1 



Where E is the expectation w.r.t. density p (x_ , ..., x ) » 



>n+l : 



'n+1 



= ( P Wl e<Xn+1 ' 9)P ll e - ' e)p e (9)de • 

The integral sign I is used in this paper to mean either integration or 
summation whichever is appropriate. Now I can be rewritten as: 



I = E 



1 



\.,ii"w - v.n"w) |*» s 
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L \ +1 li £( W> " *<i>}' 



" h v "c n+1 le »W + \ 



L E « n+1 l£ £( w> " * (i) 



We can see that 1 is minimized over all possible functions 4 { _£) when 

cb(t) = E_ t.flC ,, Y . Then I 1 = and I = I , the standard error. 
C n+1 H n+1 1 

We will denote E,. , £{£ ,.) = $(£) . 
C n+1 II n+1 

$(_g) can be, at least in theory, calculated exactly by use of Bayes rule: 



J~f(y) /p^ | (y | 9)P 1 | e (x | 0)p e (0)dedy 

$(x) = ■ 



f\\ 6 l £ I 9)p (9)de 



As it is usually difficult to employ this formula, we try to find an approximation 
to *(£) . 



t 
Here and later in this paper we use alternately <&(£) or $(x) ; the correct 

notation of course is the latter, but there is no danger of confusion. 



;| I. APPROXIMATIONS USING ORTHONORMAL BASIS 

Let X be the set of values _§_ can take. The set of Lebesque square 
integrable functions on X form a Hilbert space with measure dk (x) = p (x)dx 

f . 

I (where dx = dx 1 5 . . , , dx ) and inner product defined as 



<f 1 (l) 5 f 2 (0> = Jf 1 (x)f 2 (x)dk n (x) = Jf 1 (x)f 2 (x)p(x)dx . 

A A 

Every Hilbert space has a complete orthonormal set of functions which is 
called a Basis. 

i is °° 

Let the sequence |u.(x)i form a basis for this space then by Theory 

3 } j =0 

we know that we can write, for any function $(x) belonging to this space; 

I 

*(x) = 1 c.u.(x) where c. = <$(x),u,(x)> . 
j_Q 3 3 3 3 

? 2 

Also let <f>(x) - I c.u.(x) then the c.'s that minimize E [$(£) - <j»(_£)] 

j=0 3 3:2 1 

are given by c. = <$(x) ,u. (x)> , i.e., they are the same as above and do not 
depend on N . 

Before we go on we want to denote h(8) = E | Q f(C .,) and note that: 



by independence of 5 .-,»£ given . 

nrri 



Proposition 1 : 

Under the above notation 



c. = E e [h(e)-E c | Q u..(£)] where h(0) = E ? \ Q f ^ n +i>' 



'n+1 1 






Proof: 



- "^"' ~-4 E W^ f ^ 41 ; 



c. = <*(0,u.(^)> = E 



■m nM )*uAg 



E„E, 



e a l|6-r Ji.a^Wj 



fC^J-u^d) 



'n+l 



E 1. applies to each of the terms in the product separately since 
II 9 . 

\ ik.^W = \ + il ef(?n+l) and ^ +v ~ are independent: given 



. Thus 



(1) 



Discussion: 



c j = E e 



[\^^^'\l^m 



^iHe)'E^ Q n.(gi 



Q.E.D, 



This result is quite general since £(g ) can be "any" function (when 

n+X 



the appropriate expectations have meaning) of E, 



n+l 



K. 

Some examples are (a) f(£ ') = £ ., if we are interested in forecasting 

n+X Ttvl. 

the next claim (k = 1) or a moment of it (k > 1) . (b) f(£ . , ) = 1{| , - y) 

tttl n+X 

where 



I( w - y> 



1 W i y 



. otherwise 



To expand on (b) we note that if f(£ , ) = I(| , - y) then E i I(C , - y) 

n+l'--*- 

i'n 5 n+1 < y | CJ + 0*PU n+1 > y | II = m n+1 < y I II - ^f ^ £ y I 1,1 J = 

E e|/ [? n +i- y I e] = E ek h(0) where h(e)= p[c n+i - y I 8] ' So we are 

approximating the distribution of | t , I t at a given point y . (c) For the 

n+l -**■ 

discrete case one can set £(£ , ) = 5 (5 - y) where 



' V: ■ V- ■;■■.' ■; 



^e n+1 ~y> 



1 W 



otherwise 



then similarly as in (b) we have E | r f(5 ±1 3 = p [£ ,. = y | C] and 



*n+l' 



'n+1 



h(e) = P[C n+1 = 7 I 9] • 



The choice of a basis is not an easy task and is more or less an art in 
itself. Here we try to limit and discuss the several possibilities. 

Let a sufficient statistic for 6 exist. We call it T = T(g_) . 



Proposition 2 : 
(2) 



B Wi"W " \ +1 lT<o £( W ■ 



Proof: 



First we note that if a sufficient statistic T (E) for 9 (T and 

n — n 

a vectors) exists we can write p^i„(x j 0) = P m i„(t I 0)*H (x) where 

E 8 — ' Tun' n — 

— n' 

p„ j (t j 0) is a density function of T (£) and H (x) a function of 

ion n — n — 



can 



x r 



x alone. Thus, 

n * 



D) h +lS C|6 (X n + l^ i 0) ■ p £ +1 |0 (X n + l I 8 >'^[»fe 1 9) 



n+1 1 

P T 1 (C 4 . 1 )ie (t l (X n+l ) I 9)H l (X n+l )p T fc>ja Ct «® I *>**& 

1 n+1 n — ' 



Now: 



\+lH^ +1 



(x „+l I *) = 



w^ (Vr 



V P E ( - I Q )] 



and by (3) 



P T|8 (t l (X n + l I 6) * P T |e (t n ( ^ I 6) 
1 ' n 



^• 1} " H l (x n+1 } " E_p_, i.Ct (x) ! 6) 



n 



and 



E e p 5+i | (X n + l I 9) " P T (9 (t n« • 9) ' 
<*'« p £ ,JT (X n+l I V&" = ~" ~ eV lQ (t (xjTI) = 



^n-fl 1 *o ""* " 6'~T 10 v ~n' 

a' 



but p i Q (x I 0) = p i Q C t (x ) I 8)«H (x ) and thus (4.1) equals 

(4.2) 

Q.E.D. 

Because of the above proposition, there exists *(•) such that *(0 = ®(t(£)) . 
That is, <$(_£) is a function of the sufficient statistics for only. 

This limits our search for a basis to a set of independent functions of t , 
i.e., u . ( £) = w.(t (£_)) ■ It is our choice to write either: 

Ejw. (t(O)-w. (t(£))3 =0 density p_(x) 

£ j i — £ - 

or 

E [v. (t)*w, (t)] =0 density p (t) 
i j i i 

whichever is more helpful. 

When T is not a vector , one choice of basis is to use polynomials of the 
sufficient statistics T . This is of course possible only if the square 
integrality condition [J t'p (t)dt < « „ k integer! la satisfied. This is true, 

for example, when T e [0,°°) and p (t) is of exponential form. 

If we want to use polynomials in T as a basis, we can use a three term 
reoccurrence relation to generate them. This can be found in books on orthogonal 



'• 



expansions [4] : 



w (T) = 1 



Wl (T) = T - E T T 



w k+ i (T) = T V T) ~ a k-M w k (T) ~ WVi W) > k>1 



where 



J 



tw k (t)p T (t)dt 



<Vl 



/• 



^(t)p T (t)dt 



tw k (t)w k __ 1 (t)p T (t)dt 



k+1 



J w Li 



(t)p T ct)dt 



A similar three term reoccurrence can be shown to exist between the c. ! s by 
direct substitution of the three term reoccurrence of basis into the formula of 
Proposition 2 



E {[E , (T-w. 1 (T))]-h(e)} , , 1" 

Cj = ± -%^ * jVl - 3.c._ 2 where )| Wj (T)|| fi^ff)] 

where a. , 0, were found in calculating the basis. 

We want to note that as the calculation of a polynomial basis requires only 
the knowledge of the distribution of the collective p (x) , the calculation of 
c 's require both the distribution of 9 and T | 9 (or equivalently £ | 9) . 



II. THE ERROR IN OUR APPROXIMATION 



We said in the beginning that I = I„ + I . The best we can get (using Bayes 
rule) is I ± = . Now ^ = E ? f*(S) - <f>(I)J 2 . But 



»(i) = I c 4 » .(c) 

j=0 3 3 



and 



Thus, 



1=0 3 J 



T = P 



j=N+l J J 



Proposition 3 : 



(5) 



I 

j=N+l 



Proof: 



We kno 



w that If (£) = <*(£) .4>(0> = | c. . Also by Pythagorean Theorem 



<«(£), *u» = I c + « 

j=0 J * 



*<£>- I c u (£) 

j=0 J J 



so E (<p (i ) _ ( ^(_ i) ) 2 = I c 2 

j=N+l J 



Q.E.D. 



We want to note here that since E $ (£) = ) c. = constant, to minimize 

a. j= 3 

2 N 2 

I = Y c . is equivalent to maximize I as | c. . 

1 j=N+l J 3=0 3 

This can provide us with a way of comparison of two bases s i.e., it can 

resolve a question like which of the two bases {u.} or {v.} , that we are 

i J J 

given, is better for an N-term approximation. We will prefer the one that 
maximizes I.. . 
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TTI^ TWO TERM APPROXIMATION 

Here we restrict ourselves to two-term approximations of the form 



(I) = c U (i) + C l U l ( i ) ■ 
We pick U-Cj) = 1 and 



uAO 



g(C) - E g(£) 



V^r^ g(g) 



where g(£) any function of J_ that belongs to our space (i.e., 



Jg(x) p^(x)dx = Eg (x) < ro | . 



First, we note that u , u are orthogonal: 



E_[g(e> - Eg©] 

I _ I _ 



<u ,o > 3 - — - 

War r g(C) 



u is also normalized: 



VE^(g(U - E ? g(0) 2 

|u (0|| =- - ;== = 1 

^Var^ g(£) 



Now, using Proposition 1, 



c Q = E e h(e) 



E Q [h(0)-(E , Q g(0 - E £ g(i))] 

C l = E 9 CM8)»E |e u l(i )] = ~ " = 

— ' VVar g(u 

E Q h(e)-E , e g(0 - E e h(8).E 9 E , 6 g(C) Cov Q [h(8) ;E ? , Q g(£) ] 



\/var g(5) ^Var r g(5) 



II 



Thus, 

Cov Q [h(8);E , g(C)] 
(6) HO " E e h(9) + vag g<i1 ISW - Egg(i): 



We call 



Cov e [h(6);E e | Q g(£)] 



Thus, 



C l 



7var g(C) 

; . '"-:' From the preceding analysis, a generalization of Ericson's Theorem follows: 
(we keep the same notation) 

Theorem : 

: : " Let E, = £-, , .... £ be independent identically distributed random variables 
— 1 n 

conditional on . Let h(0) be a function of and g(£) an unbiased 
estimate of h(0) , i.e., E . g(0 = h(0) . 

If we want to approximate E i h(8) by b + ag(0 , the least square 
coefficients a , b are given by: 

Var fl h(8) 
(7) . a = -jJU^ , b - (1 - a)E e h(9) . 

Proof: 

Follows directly by (6) by setting E t g(£) = h(9) 

Q.E.D. 

Now we note that 



±2 



(8) Va ^ g ( -I> = E T Var c ) T ^W + V** T E., T g(i) 

and that if T = T(£) is the sufficient statistic of 6 for p ( x ) 6) then 
LCt(O) 8 i T g(jL) i s still an unbiased estimator of h(6) because 

E I E i T g(i) = E T | Q E ? |T,6 g( ^ ) SinCS E c|T s g( 5- ) iS a function of T only, if T 
is a sufficient statistic for 6 . Thus, E r| e E H T g(i) ;= E rlgg(S.) = &■<&■) . Also 
Var L(t(J_)) <_ Var g(£) by (8), i.e., L(t(0) has minimum variance over all 
unbiased estimators. This is an improvement over g(£) since 

2 (Var Q h(6)) 2 
C l = Var ? g(0 

2 2 
is maximized when g(£) = L(t(§)) . Thus, I = c + c is maximized. So in 

theory at least, we can find the best two-term expansion within the class of 

unbiased estimators of h(0) » 

The question is how to find L(t(|_)) .* 

1. We find an unbiased estimator for h(8) . This is readily available. 

Take f(£,) . Remember we defined h(6) = E. Q £(£) . Thus, f(£ n ) is 
i t, I b 1 

an unbiased estimator for h(0) . 

2. Find a sufficient statistic T(£) for 8 and the density function of 
T , P T | Q (t | 9) . 

3. Calculate E^f^) using P T | Q (t | 6) , p e (8) if available. 

It is obvious that Steps 2 and 3 can be very difficult in practice. Finally, we 
remark that for the exponential family the sufficient statistic is complete; thus s 
the min variance unbiased estimator is unique. 

Example I : 

We want to calculate the net premium. Then f(| , n ) = | , -, • We will use 

n+1 n-t-1 

n 

• , - i, Cl 

xcson's Theorem. Let us choose g (_§_) = 5 = which is unbiased, i.e., 
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\\e^ = E jL|e ? ' = E c|e 5 = h(e) * Then 



Var Q h(0) E [(h(6)) 2 ] - (E f) 2 

. a - - - — — -- - 



which is exactly the formula found in Blihlman [1] = 

Example II : 

.;■ For the variance principle in premium calculation, we need to approximate the 
quantity Var ? \^ n+1 



'n+1' 



Var 

^n+1 



nw - x +1 \M+J - (\ +1 uw} : 



2 
Noting that £ ,, is a function of £ , f (£ ) , we can apply previous theory 
n+1 n+1 n+i 

2 
on orthonormai basis to approximate E i | , or use a two-term expansion. 

5 n+1 Hn+l 

As for the second term, an approximation can be found by squaring the 
approximation of the net premium. 

Example III : 

We want to approximate P[£ , <_ y | _§_] . Then f (5 +1 ) = I(C + , - y) and 
h(9) = P[£ <_ y I e] = P[y | 8] . Applying Ericson's Theorem, we choose 



(D - | ! i« 4 - y) 



i-i 



which is unbiased, i.e., E , g(£) = h(6) . Then 



Var h(6) Var [P(y | 6)] 

■;a" = — B 9 " 



V "I ■«> I( „ i Il(5± - y) + Vare E i|9 i J ,«, - « 
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The denominator can be written as 



P(y | 8)(1 - P(y | 0)) 






where 






+ Var n P(y i 8) = 



EqP [y I 0] + Var Q P(y | 9) 



I P(y) - i Var Q P(y | 0) - ^ (E Q P(y | 0)) 2 + Var Q P(y | 
i P(y)(l - P(y)) + M - j\ Var P (y | 0) 



P(y) = E Q P(y J 9) 



This is exactly the formula for the credible distribution found by Jewell [7] 
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IV. A G ENERALIZATION TO MANY DIMENSIONS 



We can easily replace g by a vector g and £ by a matrix = and 



n+1 



£( W by l( ^ 3 



'n+1' 



*U> 



VW 



VW 



The relevant density is now P_i Q (X | 8) = II p_^ (x | 6) 



t=i , ye 



^n+V = n+1 W** 



m 



i=l 



£ f i ( W - h im ' 



So the problem of finding $.(=) , i = 1, . .., m to minimize I can be split to 



ra subproblems: Find cf> . ( = ) to minimize I. 



h ' \ +r , (f i ( W> "M 3 " ' i = 1 ' •••'» 



Now I can be written as I. = I,_ + I where we define 



ho ' hfc ,_ £ J (l „ + l' " *i (s > 



«„ +1 l s 



and 



I U - E_[*.(E) - ♦ <S)}' 



^n+1 ' " 
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Similar results should follow as before with the appropriate complication. 

Exampl e: 

Suppose f(| ,) = f(l n+1 ) = Kt n+1 - y) where 



-> -s- 

,1 E, ,, < y 

T /£ •*, ! n+1 — 

I(? n+1 ~ y) = , 

,0 otherwise. 



-> 

We choose to make a two-term approximation to E f (£ , ) 

~r I = n 
^n-KL ' = 



We wish to use Ericson's Theorem so we restrict ourselves to unbiased estimator 
n+1 n+1 



of h(0) = E ICC , -, ) ~ p [5 n < j I 8} » One such is readily available. Take 



? n+l |6 



1 n. 

g(H) -- J I(£. -y ) 



i=l 



Then 



a = 



Var Q h(8) Var Fff^ < y | (3; 



Var^. g{3) n 

I iU 2 " y i } i-1 ! V 

g var , n — + Var E , , ~ 

9 £ n 6 = e n 



by similar calculation as in Example III 



-j- 



Var P[y | 6] 



~ P(y)(l - P(y)) + (l - ™) Var Q P (y | Q] 
where P (^ = p[ | ^ and p[ + | Q] . pf | J , Q] § 



17 

,rAW_AC ADEMIC APPLICATION 

This paper 8 s purpose is not to advocate the use of orthogonal functions to 

i 

approximate the Bayes rule as it appears on Page 2. 

We used this approach to unite the idea of credibility in insurance whether 
this is for the approximation of the mean or the credible distribution or any other 
function as it was explained before. . 

In view of the calculations involved in approximating the Bayes rule by 
orthonormal functions, the question still exists — what are the advantages, if any, 
from the computational point of view, to calculating the Bayes formula directly. 

We can say the following: 

(a) Directly from the Bayes formula, we have to do a number of integrations 
to find i> (x) for a specific value of x_ . An improvement to that is 

to find *(t ) (a function of the sufficient statistics T ) 

n n 



$(t n } E E e|T E c |e f(5) • 

1 n ' 

Again, our effort will result in finding (with the use of a computer) 

the value of $(t ) at a specific value of t 
n n 

(h) Using orthogonal functions of the sufficient statistics for 6 , 

|w. (t )\ , we make initially a bigger effort to find them numerically. 

\ 3 n /j=0 

After that, given any function f(£ , n ) for which h(6) E E i A £(5) , we 

n+i t, j B 

use the same {w.(t )} to find the coefficients {c.} according to 
3 n 3 

formula (1) for each h(8) , our goal being to form the sum 

N 

I c.w.(t ) which is our approximation to the Bayes formula. But this 
j= 3 3 n 

sum is a function of t . So no extra integrations are needed for 

n ° 

another customer with different t . (Note n must remain the same.) 

n 
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IlTe_Exam£le 




We choose 



p |e (x I e) — 

^ ' /2tt v 



exp I- 



&J^£ =N(0,v 2 ) 



2v~ 



p Q (0) = N(y,a ) 



A sufficient statistic for 6 is T = J % . and it is easily shown that 

n i=l x 



It follows that 



P T | (t n I 6) = N(n6,nv ) 



2 2 
p (t ) = N(nu,n(no + v )) 

T n 
n 



Numerically, the first seven otthonormal polynomials in t w.r.t. density 

2 2 
P m (t ) were found with n = 10 , jj = 2 , a = 2 , v = 5 : 
n 

T! 



j 



t 


1 

t 


2 

t 








t 3 


4 
t 


t 5 


6 
t 





1. 




















1 


-1.27 


.064 


















2 


.473 


-.12 


.003 


















' .746 


.071 


.008 






.12 


-3 
x 10 








4 


-1.04 


.085 


.007 






.43 


-3 
x 10 


.53 x 1Q~ 5 






5 


.170 


-.191 


.57 x 


10" 


-2 


.52 


-3 
x 10 


.23 x io" 4 


.22 x 10" & 




8 


.854 


.082 


-.019 






.22 


x 10" 3 


.350 x 10~ 4 


.11 x 10*" 5 


.9.1 x 10~ 8 



TABLE 1: THE ORTHONORMAL POLYNOMIALS 



i^t- 
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Next we calculate the coefficients |c.l for 

1 J fj=0 

■ 
(a) hl (6) = E^^ 

(b ) h 2 (e) s e 5 | 6 c 2 

(c) h 3 (0) s P[£ < x | 9] , x = 2 . 

They are as follows : 



three different functions: 





c o 


C l 


C 2 


C 3 


C 4 


c 5 


C 6 


6 9 

j=0 J 


h 1 (6) 


1.99 


1.05 


-.156 


-.016 


-.115 


-.073 


.099 


5.135 


h 2 (o) 


10.892 


3.96 


1.43 


-.571 


-.218 


-.162 


1.596 


139.293 


h 3 (e) 


.490 


-.232 


.029 


.014 


.011 


-.051 


.065 


.3026 



TABLE 2 

Looking at the sum of squares of the c.'s , we see that in all cases the main 
< otribution comes from the first two terms. This can be used also as an indica- 

: on that a two term approximation with g(x) = t(x) would be appropriate (see 
' :tion III), Specifically, 



(a) For h 1 (9) =8 we know that the correct result is 




i=l 



10 



c Q = E e h(6) = 2 



Var h(9) 



a 2 2 x /lO 


c - 


„ , 


v n 

c. = , j > 1 . 

J 


/l ,2 2. /20 + 5 
-if n(ncr + v ) 

/a 

1 



« 1.24 



These results are very close to our estimations of c » c in 
Table 2. The rest of the c.'s are one order of magnitude smaller, 

thus their contribution to the sum of squares is minimal. 

2 2 
(b) For h„(0) = v +0 the correct result is 

9 7 22 / 2 , 2 \ 2 

it r 2 ~ 2 a. ° V a. E tlj 

n+1 n v + no \v + na / 

2 2 
setting v = 5 , a =2,p=2 



40 . 4 2 

5 - 5 + 625 C ^ 625 C 



which is also very close to what we find by using orthonormal polynomials 
(c) Finally, for h (8) =P[f; <x | 6] , x = 2 s one can show that 



\ +1 



9 o 

| T (x | t) = N(y f 9 u + v ) where 



if-^V--= .nil 

v + na' 



p 2 2 " ' 25 



and 



2 2 

2 a v 27 

11 = = — - — 

2 2 5 
v + na 



Using tables for Normal distribution, we find what appears in Table 3 under 
column "correct value". 

Now using orthonormal polynomials in t , we have: 



"' " ' ~ J L r _gL Ap P r o x i ma t i o n : 



P[£ < 2 | T = t] .78464 - „01484t 



Three-Toriri ^pP.h'^^i'il^l}. 1 



P[£ <_ 2 I T = t] .798357 - .018328t + .000087t' 

S 

Che , dculations are summarized in the following table: 



ho 


Two-Term Approximation 


Three- 


-Term Approximation 


Correct Value 


10 


.636 




.624 


.637 


20 


.488 




.467 


.500 


30 


.339 




.327 


• ,363 


40 


.191 




.204 


.242 


50 


.04224 




.099 


.1482 



TABLE 3: P[£ <_ 2 | T__ = t-J 



10 10 J 



that for big values of t we need more terms to have an acceptable 
P • iciation. 

Finally, we want to remark that the computer effort involved (not including 

• ' ng effort) was minimal in this particular example since we chose 

2 2 

: - = ° = 2 . If were big, a larger range of integration over 6 would 1 

thus increasing the computer effort vastly since we are dealing with doubi* 
\, als over x and 6 . 
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